Robust sub-shot-noise measurement via Rabi-Josephson oscillations in bimodal 

Bose-Einstein condensates 
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Mach-Zehnder atom interferometry requires hold-time phase-squeezing to attain readout accu- 
racy below the standard quantum limit. This increases its sensitivity to phase-diffusion, restoring 
shot- noise scaling of the optimal signal-to-noise ratio, So, in the presence of interactions. The con- 
tradiction between the preparations required for readout accuracy and robustness to interactions, 
is removed by monitoring Rabi-Josephson oscillations instead of relative-phase oscillations during 
signal acquisition. Optimizing So with a Gaussian squeezed input, we find that hold-time num- 
ber squeezing satisfies both demands and that sub-shot-noise scaling is retained even for strong 
interactions. 

PACS numbers: 03.75.-b, 03.75.Lm, 03.75.Dg, 42.50.Xa 
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I. INTRODUCTION 

Bimodal Bose-Einstein condensates (BECs), realized 
via double-well confinement [J-Q or internal spin states 
[^-0], hold great promise for high-precision measure- 
ments. Their long phase coherence times and controlla- 
bility of coupling and interaction parameters, make them 
ideal for the construction of atom interferometers which 
will make use of quantum correlations to reach unprece- 
dented accuracy. 

The current paradigm for matter- wave interferometers 
is the separated-pulses Mach-Zhender (MZ) scheme. The 
bimodal input state is mixed by a 50/50 beam-splitter, 
then held for a fixed duration of relative-phase acquisition 
in the presence of inter-mode bias, and then mixed again 
by a second 50/50 beamsplitter. The output number- 
difference then reflects the accumulated phase differential 
and through it, the bias present during the hold-time. 

For a two-mode coherent input, the interferometer 
phase uncertainty, A0 is limited by the Standard Quan- 
tum Limit (SQL), A9 = 1/VN, with N being the to- 
tal number of particles used. The SQL is also known as 
the 'shot- noise limit' because it is essentially the classical 
expectation for a 50/50 ensemble subject to measure- 
ments. By contrast, the preparation of strongly corre- 
lated, number-sqeezed input states, can push the phase- 
estimation uncertainty further down towards the Heisen- 
berg limit A9 = l/N, which is a factor -s/iV below shot- 
noise [1, Q , as demonstrated in two recent experiments 
0, • In these experiments controlled interactions were 
used to dynamically generate squeezing via a one-axis 
twist strategy jTo|. 

While strong interaction is essential for initial number- 
squeezing, it also limits the precision of the interferom- 
eter due to p hase-diffusion during the phase acquisition 
time [nl, [l2||. To the lowest order, phase-diffusion is the 
shearing of the initial state due to the relative-number de- 
pendent mean- field shift. Hence its rate is proportional 
to the number- variance [ll| and the phase-squeezing re- 
quired to attain sub-shot-noise accuracy, increases the 



sensitivity to interactions during the hold time. By 
contrast, states which are number-squeezed during the 
phase-acquisition period, are far more robust, but suffer 
from inherently large readout uncertainty. 

One approach to tackle this interplay between readout 
uncertainty and robustness against phase-diffusion, is to 
search for the optimal squeezing and phase-acquisition 
time which will yield the best accuracy for a given inter- 
action str eng th and particle number [12h14| | . In previous 
work [isl. Il4|. we have carried out such optimization for 
Gaussian Squeezed States (GSS), which constitute an ex- 
cellent approximation to the ground-state of the bimodal 
BEG at T=0, with the squeezing controlled by adiabatic 
variation of the interaction-to-tunneling ratio [Tsj . The 
GSS also closely approximates the states formed dynam- 
ically by the nonlinear beam-splitting of Refs. @, The 
non-interacting case was considered in [l^ and the effect 
of interactions was later studied in [ij] . In the latter work 
we adopted the view that the bimodal BEG MZ in- 
terferometer measures the 'bias', or energy differential, 
e, between the two wells, rather than the accumulated 
phase-shift 9 = eT. Thus the hold time T was taken 
as a free parameter, which can be used in supplement 
to the initial number-difference variance a to optimize 
the signal-to- noise ratio (SNR) s = Ae/e of the interfer- 
ometer, at any fixed value of N and the dimensionless 
interaction parameter 77 = UN/{2e), where U is the in- 
teraction strength. 

Using this approach, we found that the optimal prepa- 
ration changes from a number-squeezed state (trans- 
formed to hold-time phase-squeezed state by the beam- 
splitter) with ctq ~ N-^^^ in the absence of interactions 
p3| to a phase-squeezed state with (Tq ~ (here and 

below 'o' subscripts indicate optimal values) for strong in- 
teractions [ij]. The resulting optimal SNR is subject to a 
similar transition from the better than shotnoise scaling 
So ~ iV^/'^ without interactions to the strong-interaction 
behavior of So ~ \/ N/rj. Thus sub-shotnoise precision is 
lost in the MZ scheme due to interactions and the op- 
timal SNR is worse by a factor y/rj than the standard 
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FIG. 1: (Color online) Signal acquisition stage for (a) Mach- 
Zehnder interferometer and (b) the proposed Rabi-Josephson 
scheme. Husimi distributions are plotted immediately after 
the 'beam-splitting' and at the end of the hold time. Solid 
lines depict the classical, mean-field trajectories. 



magnitude of the perpendicular field and could be 
used to measure it. Such a device would thus be sen- 
sitive to parallel field gradients through variations of £ 
measured via relative-phase oscillations and to perpen- 
dicular field gradients through variations of K measured 
via relative-amplitude Rabi-Josephson oscillations. 



II. MEASUREMENT SCHEME AND 
OPTIMIZATION PROCEDURE 

As in previous work, wc use the bimodal approxi- 
mation of the two-site Bosc-Hubbard Hamiltonian [l6| . 
which accounts for the pertinent physics in current ex- 
perimental setups [l3|, 



quantum limit. 

Here we consider a different strategy in order to alto- 
gether avoid the contradiction between the requirements 
for read out accuracy (namely hold-time phase squeezing) 
and robustness to phase-diffusion (i.e. hold-time number- 
squeezing). To achieve this goal we replace the roles of 
the bias e and hopping K during the signal acquisition 
stage of the interferometer. Thus, instead of using phase- 
oscillations to monitor e we propose to directly measure 
the nonlinear Rabi-Josephson amplitude oscillation in or- 
der to determine K (sec Fig.[T]). This way, we expect that 
the preferred preparation for readout accuracy involves 
hold-time nwmfeer-squeezing instead of phase-squeezing. 
Consequently the optimal input states should not lead to 
faster phase-diffusion but in fact be more robust against 
interactions. Repeating the optimization process, we find 
that the optimal squeezing scales as CTq ^ {toNY^^ where 
To = KTo ~ and u = UN/{2K) is the perti- 

nent interaction parameter. The optimal acquisition an- 
gle To is only residually dependent on N. Consequently, 
the optimal SNR s = K/AK scales as So ~ {toN)^/^. 
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Therefore, stronger phase-squeezing ctq ^ 
ferred for stronger interactions and sub-shotnoise scaling 
So ^ u'^/^N^^"^ is maintained. 

Experimentally, current double-well setups rely on the 
combination of a 3D harmonic trap and a ID optical 
lattice, 



V{x,y,z) 



J2 "^ii^ 

q=x,v,z 



(1) 



where lOx.y.z are the trap frequencies along the respective 
axes and is the optical lattice wavenumber. Mach- 
Zendhcr intcrferomctry is implemented to measure the 
tilt Vt{z) — EzZ introduced by a field applied along the 
lattice axis. If instead a field is applied perpendicular 
to the axis, its effect would be shift the harmonic trap- 
center along X. Provided that Vb is made a;-dependent 
(e.g. by using a 2D lattice whose wavelength along x 
is long compared to the harmonic trap size or by imple- 
menting a perpendicular Vq gradient so that Vb — >■ Vq{x)), 
the barrier height and width would be sensitive to the 
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Here if, e, and U are coupling, bias, and interaction 
energies, where U > Q corresponds to repulsive inter- 
actions and vice versa. The bias e may be positive or 
negative, depending on the energy detuning between the 
two modes. The SU(2) generators Jx = {a\a2 +a\ai)/2, 
Jy = {a\a2 — d\ai)/{2i), and Jz = («i — "■2)/2, are de- 
fined in terms of the boson on-site annihilation and cre- 
ation operators , a\ , with the conserved total particle 
number ni + n2 = N = 2j . 

The standard MZ scheme consists of a fast 7r/2 beam- 
splitter rotation about J^. , followed by relative-phase ac- 
quisition during a hold-time T due to the bias detuning 
e, and an opposite tt/2 readout rotation about Jx- The 
final population imbalance j/ is used to read the accu- 
mulated phase 9 = eT from which the bias, e, is readily 
obtained. Assuming that the beam-splitter and read-out 
rotations are instantaneous with respect to the charac- 
teristic phase-diffusion time, the MZ interferometer can 
be described by the propagator 



UMZi{d,ihi) = e 



-if Jxg-ieij(i-(')/j)iz)g-if jx 



(3) 



where rj = Uj/e. The simultaneous operation of phase- 
diffusion and phase acquisition during the hold time T 
degrades the accuracy of the interferometer, as described 
inlii,lll- 

Instead, wc propose to employ a Rabi-Josephson (RJ) 
scheme, summed by the propagator 



URjiT,u,j) = e' 



r(4- 



(4) 



where u = Uj/K and r = KT, in order to determine the 
value of the coupling strength K from the final popula- 
tion difference. As in our analysis of the MZ interferome- 
ter [l^, [3] , this propagator acts on a Gaussian squeezed 
state, of the form. 
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where a is the initial number-difference uncertainty, and 
= Em=-i'3xp(-TO^/(2cr^)) w ^p2^a. Such states 
approximate well the ground state of Hamiltonian ([2]) 
with [/ > and e « [l^, as well as the dynamically 
squeezed preparations observed in recent experiments [y, 

In contrast to the MZ interferometer, where the com- 
mutation of the operators for phase-acquisition and 
phase-diffusion allows for separation of the two processes 
and the derivation of analytic expressions, we now have 
to carry out the simultaneous evolution under and 
J^. This amounts to the nonlinear mean- field Bose- 
Josephson oscillation, accompanied by the deformation of 
the initial distribution due to the variation of initial con- 
ditions. We evaluate the final population difference j/ = 
{Jz)t and its variance {IS.JIY = (Jf )t — {Jz)!-- Using 
the error propagation estimate AisT = l^Jl / [dJl / dK)^ 
we find the the SNR ratio is given by. 
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This quantity is optimized with respect to the param- 
eters, {cr, t}, in order to obtain the maximum signal-to- 
noise ratio. Because the final state |(t, t) = UbjW) de- 
pends only on the parameters {r, u, j, cr}, it follows that 
the optimal values, (cto, Tq) that give the maximum SNR, 
So = s((Jo, To), as well as So itself, are functions of the par- 
ticle number, N = 2j and the interaction-to-tunneling 
ratio u, only. Our goal is to determine the scaling of So 
with j and u and to establish if sub-shotnoise scaling of 
So is maintained in the presence of interactions. 




FIG. 2: (Color online) Optimal initial number variance CTo re- 
quired to maximize the SNR: (a) as a function of interaction- 
strength 77 at fixed j = 64, (o), 128 (□), 256 (A), and 512 
(o) for a MZ interferometer; (b) same for the proposed RJ 
scheme, with u being the relevant hold-time interaction pa- 
rameter; (c) as a function of the atom number j at fixed 
interaction strength In(Ty) = —2.5 (*), and -1-2.5 (x) 
for the MZ interferometer; (d) same for RJ scheme with same 
values of \n{u) using the same notation convention. Sym- 
bols in all panels denote full numerical simulation results. 
Lines in (a,c) correspond to analytic expressions obtained 
in Ref. [lj| for weak interaction (dashed) and strong inter- 
action (solid). Solid lines in (b,d) correspond to the BBR 
truncated-cumulant-expansion [ig |. Dotted lines in panels 
(c),(d) mark the interaction-free optimal-squeezing scaling 
(To = j^^^ whereas dash-dotted lines mark the coherent-state 
(classical) variance a — ^/j 



III. NUMERICAL OPTIMIZATION RESULTS 

In Fig. [2] we plot the results of a numerical evaluation 
(solid lines) of the optimal squeezing ctq required to ob- 
tain the best SNR in both the MZ (Fig. ^,c) and RJ 
(Fig. [l]3,d) schemes, as a function of the pertinent inter- 
action parameters {1] for MZ or u for RJ) and particle 
number j. Symbols correspond to the analysis oulined 
below in Section IV. For both MZ and RJ schemes, the 
optimal squeezing, as well as the resulting SNR (Fig. SI, 
coincide in the weak interaction regime rj, u < 1/2, with 
the previously studied interaction-free MZ propagation 
p^ . This is expected because without interaction the 
choice of axes is immaterial and the RJ sequence may 
be viewed as a simple rotation of the standard MZ in- 
terferometer. For stronger 77, m > 1/2 interactions (corre- 
sponding to the Josephson regime [l8| for the RJ scheme) 
however, significant differences are observed. 

For the MZ sequence, the optimal initial number vari- 
ance increases with increasing interactions, until a transi- 
tion from optimal initial number-squeezing (cto < ■\/j/2) 
to optimal initial phase-squeezing (cto > ■\/j/2) takes 
place as rj crosses unity. This transition results from 
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FIG. 3: Optimal acquired angle 60 = eTo, To = KTo required 
to maximize the SNR for (a) the MZ and (b) the RJ schemes, 
respectively, as a function of the pertinent interaction param- 
eter 77, It and the particle number j. 



the said contradiction between projection noise min- 
imization (initial number squeezing, hold-time phase- 
squeezing) and phase-diffusion control (initial phase- 
squeezing, hold-time number-squeezing). Thus, while for 
weak interaction (77 ^ 1), we obtain sub-Poissonian scal- 
ing (To fx j^^^ [Tsj . in the presence of strong interactions 
(77 > 1) wc have super- Poissonian optimal- input num- 
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FIG. 4: (Color online) Optimal precision Po = InSo. Panels 
are arranged in the same order as in Fig. [2] Values of j in 
Panels (a),{b) are also the same. In (c) and (d), the values of 
\n{rj) and ln(u), respectively, are -2.5 (*), (+), 2.5 (x), and 
5 (■). Using the same notation as in Fig. (2] markers denote 
numerics, whereas lines are analytic estimates (a,c) and BBR 
approximate values (b,d). Dotted lines in (c),(d) mark shot- 



noise scaling s ' 
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ber distribution with (Tq oc y/rjN [14| . In comparison, 
the optimal squeezing in the RJ scheme does not follow 
a similar compromise. Initial number squeezing is pre- 
ferred for both weak and strong interactions, and in fact 
becomes stronger with increasing the value of the inter- 
action strength u. The optimal squeezing in this case, is 
found numerically to scale as {toNY^^ . 

Significant differences between the standard and pro- 
posed schemes also appear in the behavior of the opti- 
mal acquired angle {9o = eto for MZ and Tq = Kto for 
RJ), shown in Fig. [3l For neg ligible interactions, both 
are optimized by ~ Stt/S [IJI as the pure Rabi oscil- 
lation is a simple rotation of the phase-oscillations of 
the Mach-Zehnder interferometer. However, the j de- 
pendence of the optimal t in the RJ scheme in the pres- 
ence of interactions is much weaker than that of 9o for 
the MZ interferometer. The resulting maximized preci- 
sion po = Inio So {po corresponds to the number of sig- 
nificant read-out figures) is plotted in Fig. SI using the 
same conventions and ordering. For the MZ scheme, the 
weak-interaction sub-shot-noise scaling So oc j^/^ [1^ is 
replaced by So oc yN/rj for ?/ > 1, i.e. worse than the 
standard quantum limit. The RJ scheme is far more ro- 
bust to interaction and the SNR deterioration is slower 
than So oc u~^/^ as compared to Sq oc u~^I'^ for the MZ 
scheme [l^ . Consequently, the RJ scheme retains sub- 
shot-noise scaling, with Sq ^ {toN)^^^. 



IV. BOGOLIUBOV ANALYSIS 

For the MZ scheme we were able to derive exact ana- 
lytic expressions for the dynamics with general Gaussian 
initial conditions and extract the pertinent scaling laws 
for (Jo and Sq by retainin g lea ding order terms in the var- 
ious interaction regimes jl2l fl^ . Such derivation is not 
possible for the proposed RJ scheme due to the incom- 
mutability of the Jx and terms during the acquisition 
stage. In order to analytically prove the numerically- 
observed scaling relations of the RJ scheme, we employ a 
density-matrix cumulant-expansion technique developed 
under the name Bogoliubov Back- Reaction (BBR) pij| . 
The hierarchy of dynamical equations for the Ji oper- 
ators, is truncated at second-order to obtain equations 
of motion for the mean-field single-particle Bloch vec- 
tor R = (J) / j and the correlation functions Aij = 
{{,liJj + JjJi)—2{Ji){Jj))/j'^. Neglecting the back-action 
terms, these read, 

Rx = —2uRyRz, 



Ry - 

R. = 



^Xll 



+ 2uRxRz, 

-Ry 

-iu {RyAxz + RzAxy) , 

2Ayz + 4u (RxAyz + RzAxy) , 

+ 2uRz {Axx - Ayy) + 2u (R^Axz - RyAy,) , 
-Axy + 2u {RyAzz + RzAyz) , 
Azz- Ayy + 2u{RxAzz + RzAxz), (7) 



where / = df /dr. 

Given the values of R and A.^ after the preliminary 
tt/2 rotation about J^, 

Ro = (0, 1, 0) , Axy,o = Axz,o = \z,o = 0, 



2a2 



A 



yy,o 



16(7" 



^. (8) 

J 



we obtain an exact solution to Eqs. ([7]) in terms of ellip- 
tical Jacobi functions (see Appendix). To leading order, 
we find that, 



(AJ/)V/ = Azz(r)/2 



36cr2 32a 



-7- (9) 



(9J/ 



\Ry{T)\ « y/l-T^-u^T^. (10) 



Substituting into Eq. ([6]) and focusing on the strong in- 
teraction regime itr 3> 1, we obtain. 



s'{T,a) 



(1 - u't^ 



(11) 
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Optimizing s(r, cr) with respect to r and cr, results in 

^ A{u, r) , ^2 = B{u, t)2'I''t''I\ (12) 
where the functions 

only depend weakly on u and r and could therefore be 
assumed constant. Thus we conclude that the scaling of 
the optimal initial number variance CTq and hold time Tq 
required to obtain the best SNR s in the RJ scheme, is 



functions, 



(13) 



This is consistent with the numerically obtained ctq ~ 
y-i/6ji/3^ ^ j2/3y-i/3 scaling in Figs. EH The re- 
sults of the BBR optimization are shown as symbols in 
Figs. [2jb),(d) and Figs. |4l^b),(d) and agree well with the 
numerics. 



V. CONCLUSIONS 

Due to their strong interaction, phase-diffusion is an 
eminent obstacle to the realization of sub-shotnoise inter- 
ferometry in dilute quantum gases. The standard Mach- 
Zehnder scheme borrowed from linear optical interferom- 
etry, suffers from an inherent contradiction between the 
preparation of a number-squeezed input state required 
for sub-shotnoise precision, and its increased sensitivity 
to interactions during the signal acquisition time after 
it has been converted to a phase-squeezed state by the 
beam-splitter [l2|, [T^] . Placing the measured perturba- 
tion between the sites so it affects the odd-even detun- 
ing instead of the bias between sites, and monitoring 
Rabi-Josephson amplitude oscillations instead of phase- 
oscillations, removes this contradiction. We have shown 
that number-squeezed states satisfy both readout and ro- 
bustness requirements and that in fact increased number 
squeezing optimizes the input states in the presence of in- 
teractions. This is distinct from the MZ scheme, where a 
transition occurs from optimal hold-time phase-squeezing 
(optimizing readout) to optimal number-squeezing (op- 
timizing robustness) as the interactions cross a critical 
magnitude [ij]. As a result, the best SNR ratio obtained 
by the RJ scheme retains sub-shotnoise scaling, as op- 
posed to the super-shotnoise scaling of the SNR of a MZ 
interacting atom interferometer. 



Appendix A: Solution of the BBR equations 

Given the set of equations (0 with the initial condi- 
tions ([5]), an exact analytic solution for the mean- field 
equations for R is found in the form of Jacobian elliptic 



Rx — uRI, 

Ry = ^/^-R'i- u^Ri, 



Rz = —Rm.cn I 



(T(l+4u2)i/^„i^(fc)) , (Al) 

where i?^ = 2/(1 + Vl + 4u2), p 

(1-1/V1 + 4m2) /2, i^(fc) d0 (1 - fc2 sin^ (f,y'^^ 

is the quarter period, and 



The equations for the fluctuations A^j, may then be 
solved in the form. 





Axy 


A. 


^xy 


Ayy 


A. 


A,, 


\z 


A' 



Axxfl 











Ayyfi 











Azz,0 



(A2) 

where the elements Cij of the orthogonal matrix C satisfy, 

Cla = -2wi?^C2a ~ '^uRyCsa, 

C2a = (1 + 2uRx)c3a + 2ui?^Cia, 

C3a = -C2a, (A3) 

with initial values Cab{t = 0) = Sab and a,b = 1,2,3. 
Parametrizing the solutions in terms of the variable p ~ 
—Rz and noting that, 

dCab dCab 



dT dp 
we find the following series expansions: 



(A4) 



cii = 1 - 2uy^anP 



2n+2 



C21 



p2 _ u2p4 



5](2n+l)«„p2„+i^ 



C31 = anP' 



2n+l 



00 

^2ri-|-l'^ 



C22 = Vl-p'-wV4^/3„p'", 

00 „ 
— \ ^ ^" n2"+l 

00 

C13 = -2u^7„p2"+i , 



Tl = 



C23 = -v/l-p2-ii2^4^2n7„p2»-i, 

00 

C33 = E^"/''" • (A5) 



n=0 
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The coefficients q;„,/3„,7„ are given by tlie recursion re- 70 = 1, 71 = — 1/2 
lations, 

u 2u 2n-l 2 - 3n + 2 

ai - o ' "2 - -jT ' In+i = — -TT7n + u — — — -- — — --7„-i, n>l. 

o ^o 2n + 2 (2n+l)(n+l) 

2n 2 2?^ — 3 

Q^n+i = 7; — T^Qf" + " n — —^ctn-i, n>2, Acknowledgments 
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